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ABSTRACT 

We discuss an optimal .R-parity breaking SUSY solution to the Rb excess as 
well as the ALEPH 4-jet anomaly. The latter arises from the pair production 
of stop via chargino decay at LEP1.5, followed by its i?-violating decay into a 
light quark pair. The model satisfies top quark and Z-boson decay constraints 
along with gaugino mass unification. 
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Two of the intriguing results from LEP which have attracted a good deal of theoretical 
interest are the Rb and the ALEPH 4-jet anomalies. The first anomaly refers to the LEP1 
value of Rb(= r^/r| ad ) being ~ 2a larger than the SM prediction |T], Q|. The second refers 
to the anomalous 4-jet events recently reported || by the ALEPH experiment at LEP1.5, 
each of which seems to consist of dijet pairs with a common invariant mass ~ 55 GeV. 

It is now widely recognized |fj. [5]. that the minimal supersymmetric standard model 
(MSSM) offers a viable solution to the Rb anomaly in the low tan f3 region if one assumes a 
relatively light top squark (t\) and chargino (xi). The assumption of a -R-parity violating 
Yukawa interaction term would invalidate the canonical missing {$t) signature for 
superparticle production and thus the LEP1.5 limit on chargino mass [0], 

m~ > 65 GeV . (1) 

In this case, pair production of charginos at LEP1.5 can offer a possible explanation for the 
ALEPH 4-jet events as was recently suggested in refs. |8], |9|, |i~0| . In particular, Chankowski 
et al. [It| have discussed a variety of such i?-parity breaking SUSY solutions to the ALEPH 



anomaly, which can also account for the Rb excess. The purpose of this note is to focus on 
what appears to be an optimal i?-parity breaking SUSY solution to the Rb and ALEPH 
anomalies, within the constraints of top quark and Z-boson decays as well as that of 
gaugino mass unification. We shall see below that it can quantitatively account for the 
essential features of the ALEPH 4-jet events as well as for the Rb anomaly. 

Explicit breaking of i?-parity introduces additional Yukawa terms in the superpoten- 
tial UTTJI 

Wft = IXijkL^Et + K^UQ.Dl + \x' jk UtD]Dl , (2) 

where Q, L (U, D, E) denote the quark and lepton doublet (singlet) superfields and the 
subscripts denote the generation. Symmetry considerations imply that Ay& = —Xjik and 
\"j k = —X" k j- Proton stability demands that all products of the form (A'A") be vanishingly 
small and this, conventionally, is ensured by stipulating that either the baryon number 
violating couplings (A") or the lepton number violating couplings (A, A') are non-zero but 
not both. 

Let us briefly discuss the various i?-parity breaking SUSY scenarios that have been 



suggested as explanations of the 4-jet excess. The first of these |[L2|| , which, in fact, pre- 
dates the anomaly, assumes the lightest superparticle (LSP) to be a sneutrino instead of 
the neutralino. Pair production of sneutrinos (e + e~ — > z>*z>) followed by their diquark 
decays — through one of the A' couplings — can then lead to a 4-jet final state. An ade- 
quate production cross section can be obtained for v e (via a light, and gaugino-dominated, 
chargino exchange), but the other event characteristics have not been analysed so far. The 
second scenario |13| suggests pair production of left-handed b squarks, followed by their 
diquark decays through a A" coupling. However, a light bi in the required mass range 
(~ 55 GeV) seems to be disfavoured by the precision measurements of electroweak observ- 
ables |fL0[l. More recently, there have been two suggestions based on pair production of 
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charginos with mass ~ 55 GeV [§, [J . One of them || assumes the decay sequence 



Xl^W*x°i, W*^qq', x?^gig 2 g 3 , (3) 

where the star denotes an off-shell W boson. Although this is the most natural scenario in 
terms of the MSSM mass spectrum, it seems to be disfavoured on several counts. It shows 
much broader distributions in the difference of the dijet masses as well as their sum than 
the ALEPH data ||. Moreover, the leptonic decay of one of the W*s would imply roughly 
as many anomalous events with an isolated lepton as without it, which could not have 
been missed. Finally, one would have a significantly large fraction of the events with more 
than 4 jets if one applies the ALEPH jet algorithm. The other suggestion || assumes the 
chargino decay 

Xi rqq' (4) 

to dominate over its i?-conserving decay into the LSP. In this case, the outgoing r would 
be too hard to have been missed, unless one assumes the exchanged sneutrino mass to be 
very close to that of the chargino. It may be noted here that none of the above scenarios 
address the issue of the Rb anomaly, to which we now turn. 

We start by considering the standard /^-conserving MSSM solution to the Rb anomaly, 
within the constraints of top quark and Z boson decays 0, |J. Gaugino mass unification 
shall be assumed all along as it is very closely related to the successful MSSM prediction 
for the unification of the SU(3) ® SU(2) ® U(l) gauge couplings. Thus, the masses of the 
corresponding gauginos g, W and B are related via 

M 3 = — sin 2 9 W M 2 ~ 3.5M 2 and M x = - tan 2 8 W M 2 ~ 0.5M 2 . (5) 
a 3 

The physical gluino mass is related to the running mass through the QCD correction 
factor p]l 

nig = (l + 4.2— ■) M 3 ~ 1.15M 3 ~ 4M 2 . (6) 

Thus, a single gaugino mass (M 2 ) along with the higgsino mass parameter (/i) and the 
ratio of the two Higgs vacuum expectation values (tan/5), determine the gluino mass as 
well as the masses and compositions of the two chargino and four neutralino states [JT5f] 
i.e., 

xf L = VaW£ + V a H£, fa = U a W£ + U a H±, 

^ ~ ~ (7) 
X° = N tl B + N l2 W 3 + N i3 H? + N tA H° . 



In the scalar sector, the large Yukawa term for the top results in a mass hierarchy and 
thus the lighter (and predominantly right-handed) stop, 

ii = cos 6i i R + sin 9f II , (8) 

is expected to be significantly lighter than the other squarks. We shall be primarily inter- 
ested in this stop. 
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In the low tan/3 region of our interest, the SUSY contributions to Z — > bb arise from 
the triangle graphs involving xOifik an d UtjXk exchanges as well as the tiXj loop insertions 
in the b and b legs [[5], || The b vertices are dominated by the btfiXi Yukawa coupling 

V 2m w sin p 

which favours large V12, i.e., the higgsino-dominated region <C M2). On the other 
hand, the ZxoCj couplings 

0% = - 1 - (cos 26 w 6 ij + UaU n ) and Og = -~ (cos 29 w 6 ij + Ki^i) , (10) 

favour large Un and Vn, i.e., the gaugino-dominated region ^> M 2 ). (The corre- 
sponding Ztiti coupling is suppressed by the U(l) coupling factor sin 2 ^.) The need 
for sizeable b as well as Z couplings then implies that the largest SUSY contribution to 
Rb (5Rf VSY ) occurs for the mixed region (|/z| ~ M2) — -corresponding to a 7-dominated 
LSP — rather than for the higgsino-dominated region |5|, Moreover, it seems to favour 
negative /i over the positive \i region ||. We shall, therefore, restrict ourselves to the 
former. We shall consider 

5Rf VSY ~ 0.0020 — 0.0025 (11) 

to be a viable solution to the Rb anomaly. It would exactly account for the discrepancy 
between the current experimental value of _R^ xp = 0.2178 ± 0.0011 and the standard 
model value (for m t = 175 GeV) of Rf M = 0.2156, as well as close the gap between the 
a s (m 2 z ) estimates from LEP1 and from deep inelastic scattering. 

Low masses for t\ and xii required for a suitably large <5-R^ USY , may, however, result in 
significant new decay channels for the top quark, viz., t — > ^1X1,2- The decay amplitudes 
are dominated by the fiXi^m an d inX^iL Yukawa couplings 0, ||, [L7| : 

^^^i-cos%(sin%), (12) 
mw sin p 

with i = 1,2. The higgsino-dominated region corresponds to large components in Xi 21 
and hence a large SUSY Branching Ratio {B~s) f° r top decay. This quantity is relatively 
small in the mixed region since only x\ nas a l a rge component while Xi — 7- Thus the 
upper limit on Bs from the CDF top decay data [0 favours the mixed region as well. We 
shall take a rather lenient value for this limit || Bl : 

B s < 0.4 . (13) 



Recently, a systematic scan of the MSSM parameter space was carried out to obtain 
the best SUSY contribution to Rb within the constraint of top quark decay. The optimal 
value of 5Rf VSY is obtained at small negative value of the stop mixing angle (% ~ —15°) 
and small stop mass 

m h - 50 — 60 GeV , (14) 
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Table 1: SUSY contributions to (5Rf VSY ) and to the top BR (Bs) are shown, along 
with the light chargino and neutralino masses, for three representative points in the (M 2 , n) 
plane. (All masses are in GeV and 9$ = — 15°.) For each set, the top row corresponds to 
the allowed MSSM parameter space in the R-conserving scenario. 



which lies in between the LEP1 limit (m^ > 45 GeV) |19| and the DO excluded region 
(m^ 7^ 65-85 GeV) PD| . LEP1.5 imposes no additional bound on as the pair- 
production cross section for the right-handed stop is small. 

In Table 1 we display the phenomenological consequences for three representative 
points in the (M%, /V) plane : 

A : (150, -40) GeV, B : (60, -60) GeV and C : (40, -70) GeV, (15) 

belonging to the higgsino-dominated (A) and mixed regions (B, C). The points are chosen 
close to the LEP limit so as to give the best values of SRf VSY in the respective regions. 
For the higgsino-dominated point (A), m~ is close to its LEP1.5 limit ([]]). Still 5Rf VSY 
is smaller than the required value (|Tl"|), while the contribution to top BR exceeds the limit 
(|T3|). On the other hand, the mixed region points (B, C) are seen to give viable values of 
5Rf VSY , while satisfying the Bs limit. Note that, in this region, the chargino mass is safely 
above the LEP1.5 limit (|l|); the most important constraint is rather set by the requirement 
that 

m ~o + m ~o ~ m z (16) 
which follows from the stringent LEP1 bound 

BR(Z - xl X§) < 5 x 10~ 5 (17) 



deduced from the negative search for acoplanar jets using the ^signature of Xi The 
best values for 5Rf VSY and Bs are obtained for point (C). However, it corresponds to a 
gluino mass of 160 GeV, barely above the Tevatron lower limit PT . 



We now turn to effects of i?-parity breaking via the A" couplings. As the LSP (xT) 
now undergoes hadronic decay, the .^signature is n o longer applicable. Thus, the LEP 1.5 
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bound ([!]) is inoperative. Moreover, the LEP1 bound ( fPTD on Z decay into neutralinos is 
now replaced by 

£r(Z^4) <3MeV (18) 

corresponding to the la error in T^ ad . This is weaker than (|17]) by more than an order of 
magnitude. Nonetheless, as we shall see below, it acts as a strong constraint on efforts to 
reduce the chargino mass below 65 GeV in the mixed region. 

The goal then is to have m~ < 65 GeV, while satisfying (^) so that chargino pair 
production at LEP1.5 can be a source for the anomalous 4-jet events. In the higgsino- 
dominated region this is achieved most easily by decreasing while in the mixed region it 
can be achieved only by increasing tan (3. The former also has the advantage of increasing 
5Rf VSY simultaneously. Table 1 shows that it is possible to go down to m~ = 58 GeV 
within the constraint (|T8"D in the higgsino-dominated region. The stop mass can then be 



reduced to 50 GeV so as to allow the two body decay mode 

Xi -> hb , (19) 

which shall be assumed later on. Reducing the chargino and the stop masses has the effect 
of increasing 5Rf VSY to the respectable value of 0.0019. Unfortunately, Bs is untenably 
large. 

For point B (|yu| = M2), it is not possible to drive down vn~ x below ~ 65 GeV by either 
of the above methods, while still satisfying flT8"|). However, for the mixed region point 
C — which offers the best value for SRf VSY — it is possible to achieve this by increasing 
tan (5. As Table 1 shows, one can go down to m~ = 56 GeV, within the above constraints, 
by increasing tan/? to 2.6. Decreasing the stop mass to 50 GeV ensures the two-body 
decay (|19|) with a soft b. Note that the large value for SRf VSY is obtained within the Bs 
constraint (0). Also, the value of tan/5 is now more reasonable. Furthermore, the gluino 



mass limit from the Tevatron is no longer applicable. 

We shall quantitatively pursue the i?-parity breaking SUSY scenario summarised in 
the last row of Table 1 as a possible solution to the anomalous 4-jet events. The pair 
production of charginos, followed by their two-body decays QT9D, results in tit* along with 
a soft bb pair. For 

m ~ = 56 GeV and m h = 50 GeV , (20) 

the b momenta are always less than 6 GeV and so the 6s are expected to largely miss the 
lifetime tag. Furthermore, leptons from 6-decay do not survive isolation cuts. The stop 
can decay either directly 

tlX+ds, (21) 
or through the i?-conserving loop process 

h - cx\ , (22) 

followed by 

Xi *Luds ,cds . (23) 
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Figure 1: Pair production cross section for a 56 GeV chargino as a function of the 
sneutrino mass for the two different LEP1.5 center- of -mass energies. 



We do not consider stop decay modes with b quarks in the final state since these will lead 
to a large number of 6-tags in conflict with the ALEPH results ||. The loop decay (|22"D 
is a third-order electroweak process and hence has a very small width ||22|| . Consequently, 
the i?-violating decay (|21|) dominates over a very large range of the Yukawa coupling 

Kds ^ 5 x 10- 5 . (24) 



On the other hand, one requires 



~ 5 x 10 ( 25 ) 



for the LSP decay to occur within 1 cm |23j]. Note that the direct decay fl2"T| ) dominates 
as long as \" ds is larger than the relatively modest limit of (^) irrespective of the other 
Yukawa couplings. Thus the direct decay is at least as natural [[25[] as the alternative route 
of ( p2| & We shall see below that the former can quantitatively account for the ALEPH 
events, while the latter cannot. 

Figure [I] shows the chargino pair-production cross-section at the (LEP1.5) energies of 
130 and 136 GeV. Since the interference between the s-channel ('j/Z) and the t-channel 
(z/ e ) is a destructive one, the cross section increases strongly with the sneutrino mass. For 
the rest of the analysis we have averaged the cross sections at the two energies assuming 
a sneutrino mass of 200 GeV. 



We have studied chargino pair-production and subsequent decay via the stop (|T9| ) using 
a parton level Monte Carlo program. Both the two-body (pTJ) and the four-body ([22] & 
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3p decay modes have been considered. To estimate the effect of jet energy resolution, 
we have compared the results with and without the suppression of soft partons (arising 
mainly from 6-decays), having energy less than the ALEPH resolution error M 



a E = (0.6 J E( GeV) + 0.6) GeV (1 + cos 2 



(26) 



We found the difference to be small. The results presented below correspond to the sup- 
pression of the soft partons having energy < 
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Figure 2: Event distribution as a function of the number of jets after the initial partons 
have been clustered with Durham/ JADE algorithms $1% Both decay modes of the stop 
are shown. 

The parton jets are merged applying the Durham algorithm till 

y Dur = 2 min(£ 2 , E]) (1 - cos%)/s > 0.008 . (27) 

Events merging into less than 4 jets are reclustered with the JADE algorithm till 

7/jade = 2E i E j (1 - cos %)/s > 0.022 . (28) 

Figure |2| shows the resulting distributions in the number of jets for the two decay modes 
in question. For the direct decay mode, the 4-jet sample dominates in agreement with 
ref.[[|, while for the neutralino- mediated case the number of 5-jet events is uncomfortably 
large. In the former case, QCD radiation effects could result in a few 5-jet configurations 
as observed in 0]. 

The 5-jet events of Fig. |^ are then clustered down to 4 by merging the two jets with 
the smallest invariant mass. Within this sample, events having the smallest dijet invariant 
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mass < 25 GeV are rejected ||. This reduces the number of events by ~ 40%, which is 
compatible with ref.||. Figure |3|(a) shows the distribution of the remaining 4-jet events 
in the minimum difference of the dijet invariant masses. At this stage, both distributions 
are in reasonable agreement with ref. 0. QCD radiation effects are expected to cause a 
marginal broadening of these distributions. 
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Figure 3: Event distribution (after clustering down to 4 jets) and rejecting events with 
smallest dijet invariant mass < 25 GeV. (a) As a function of the minimum difference of 
the dijet invariant masses, (b) As a function of the sum of the dijet invariant masses for 
the pair with the smallest mass difference. 



Figure |3](fr) shows the corresponding distributions in the sum of the dijet invariant 
masses (the pairing decided by the minimum difference). For the direct stop decay (plj) , 
this distribution is sharply peaked at 105-110 GeV, in agreement with ref. [Q. The slight 
downward shift of the peak from 2m ~ is due to the suppression of the soft jets (E < cr#) as 

discussed above. On the other hand, the four-body decay of stop via \i (HI & HD is seen 
to result in a very broad distribution. Finally, it should be noted that the normalization 
of the solid curves in Figs.|3](a,6) are about twice as large as the ALEPH event size (~ 9). 
One could reduce our event rate by assuming a smaller sneutrino mass (~ 150 GeV). On 
the other hand, several of the ALEPH cuts such as the number of charged tracks and 
individual jet masses etc., could not be incorporated into our parton level Monte Carlo. 
As these will, typically, result in a loss of efficiency, we leave this excess in normalization. 

In summary, relatively light stops (mj ^ 60 GeV) and charginos offer a viable MSSM 
solution to the R b anomaly within the constraints of top quark decay and gaugino mass 
unification. Assuming an i?-parity violating Yukawa coupling A" in the superpotential, 
it is possible to bring down the chargino mass below 60 GeV as well, while respecting 



8 



the Z-decay constraint. Consequently, the pair production of stops via chargino decay at 
LEP1.5 can offer a viable solution to the ALEPH 4-jet anomaly as well. The direct decay 
of stop into a light quark pair is expected to be the dominant mode for \" ds >5x f CT 5 . 
This can account for the essential features of the ALEPH events at a quantitative level. 

We gratefully acknowledge discussions with Sunanda Banerjee, Piotr Chankowski, Atul 
Gurtu and Stefan Pokorski, and computational help from Sreerup Raychaudhuri. 
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